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Abstract. We calculated diffusion timescales for Ca, Mg, Fe in hydrogen atmosphere white dwarfs with temper- 
atures between 5000 and 25000 K. Wi th these timescales we d eter mined accretion rates for a sample of 38 DAZ 
white dwarfs from the recent studies of lZuckerman et al.l (120031) and lKoester et all i2005l) . Assuming that the ac- 
cretion rates can be calculated with the Bondi-Hoyle formula for hydrodynamic accretion, we obtained estimates 
for the interstellar matter density around the accreting objects. These densities are in good agreement with new 
data about the warm, partially ionized phase of the ISM in the solar neighborhood. 
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1. Introduction 

Gravitational separation of elements in the strong gravi- 
tational field in the outer layers of white dwarfs has been 
known for a long time to be the primary phy sical process 
deter mining the atmospheric composition l|Schatzmanl 
Il947|) . In the absence of competing processes (stellar wind, 
radiative levitation, convection), the heavy elements dif- 
fuse downward on a timescale that is always short com- 
pared to the evolutionary timescale, leaving the lightest 
element that is present floating on top. In most cases this 
is hydrogen, leading to the spectral type DA; if no hydro- 
gen is present, the atmosphere is pure He (spectral types 
DB, DC). 

However, there are some exceptions to this rule. In 
cooler white dwarfs below 25000 K, which are the topic of 
this study, we sometimes find traces of Ca (through the 
strong Call resonance lines) , occasionally accompanied by 
Mg and Fe. For He-rich atmospheres this has been known 
from the time of the first white dwarfs discovered, since the 
well known vMa 2 is a member of the class now called DZ 
(Z always stands for the presence of metals, DZ means no 
other elements visible). The class of hydrogen-rich objects 
with metals (DAZ, with visible hydrogen and metal lines) 
has really only come into existence over the past 8 years 
with the detection of a few dozen m embers in large sur- 
veys at the Keck and VLT telescopes JZuckerman ■V Heidi. 
Il998l IZuckerman et all l2003t iKoester et all 120051) . The 
reason for this difference is purely observational bias; be- 
cause of the much lower transparency of a H atmosphere 
at these temperatures, the Call equivalent widths in the 
DAZs are about a factor 1000 smaller than in the DZs at 
the same abundance. 
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The diffusion timescales for the heavy elements to dis- 
appear from the atmosphere - or in case of an outer con- 
vection zone from the bottom of this zone - is always short 
compared to the cooling timescale. A possibly competing 
process, radiative levitation, becom es completely negligi - 
ble for Ca below T cff = 25000 K l|Chaver et all 1199.4 
This implies that the observed metals cannot be primor- 
dial. They must have been supplied from the outside, and 
currently the most widely accepted mechanism for this is 
accretion of interstellar matter. There are still some prob- 
lems with this scenario, the most serious being the absence 
of co ol, dense clouds in the ISM of the solar neighbor- 
hood. lAannestad et alJ l|l993j) tried to correlate the posi- 
tions and motions of DZ (helium-rich) stars with condi- 
tions of the ISM, but with inconclusive results. Alternative 
explanations are therefore still b eing considered by severa l 
authors (see e.g. the discussion in lZuckerman et allfeOOfljh 

Because the timescales for diffusion and thus metal 
visibility in the DZs ty pically are on the order of 10 6 yr 
jPaauette et all fl986bl see also Sect.0}, the stars could 
have traveled a large distance (of order 50 pc) since the 
accretion ended, so that a lack of correlation is not too 
surprising. This situation has changed dramatically with 
the recent detection of a large number of DAZs. In these 
hydrogen-rich atmospheres the diffusion timescales are 
about 3-4 orders of magnitude shorter, practically mean- 
ing that we can only expect to observe metals if the accre- 
tion is still going on at the present time. The DAZs thus 
offer two advantages compared to the (helium-rich) DZs 

— we know where the accretion occurs and can compare 
directly with the conditions of the ISM in that region 

— we argue that we very likely observe a steady state ac- 
cretion/diffusion, meaning that we can determine the 
accretion rate responsible for the observed abundances 
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Fig. 1. Logarithmic Ca abundances (circles) and upper limits(crosses) for the combined sample (see text). The con- 
tinuous line indicates a constant equivalent width of EW=15mA, taken as observational limit for most objects 



There is thus a much more promising chance to learn 
something about the accretion process and/or the condi- 
tions in the ISM from an analysis of the observed metals 
in the DAZs. 

2. The accretion-diffusion scenario 

The idea that metals are accreted onto white dwarfs and 
then diffused downward out of sight has been developed 
by many authors. It was finally explored in great detail 
and transformed into a clear and testable theory in a se- 
ries of three fundamental pap ers by the Montreal group 
l|Dupuis et all Il992t Il993albll . The basic idea of the so- 
called "two-phase accretion model" is that the white dwarf 
travels for most of the time in the hot tenuous phase of the 
ISM with negligible accretion rates. After typically 5 x 10 7 
yr, it enters a dense cloud, where the accretion rate is high 
and the metal abundance in the outer layers approaches an 
equilibrium value between accretion and diffusion. After 
about 10 6 yr traveling through the cloud, the star emerges 
and the metal abundance decreases exponentially with the 
diffusion timescale. 

The observations available at that time - almost ex- 
clusively of DZs and DBZs - agreed quite well with 



the predictions of this model, although it obviously is 
an extreme oversimplification as emphasized by the au- 
thors. It was subsequently called into question with 
the findings that the required dense neutral clouds, 
which would be detectable through Nal absorption, 
are not found within th e Local Bubble, where these 
white dwarfs are located llSfeir et all Il99ft IWelsh et all 
I1998L Il999t fRedfield fc Linskvl 120021 iLehner et all 120031) . 
Nevertheless we will keep this model as our first hypothesis 
for this study since it is specific enough to make quanti- 
tative, testable predictions. 



While dense, cool gas seems to be absent in the Local 
Bubble, the presence of warm, partially ionized gas has 
been c learly demonstrated by iRedfield & Linsk vf ir2002l 
20043). The average number density of neutral hydrogen 
is around nni = 0.1 cm~ 3 , but very little is known about 
the spatial distribution. It is also unclear whether accre- 
tion und er these conditions wo uld follow the Bondi-Hoyle 
rate (see iKoester et all l2005|) . However, following argu- 
ments explained below, and in order to have a different 
model for testing against the standard assumptions, we 
will define our alternative hypothesis: 
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Table 1. Data for 38 objects with detected photospheric 
Ca. Atmospheric parameters T iff, log a, and logarithmi c 
Ca/H ratio, JCa/Hl , are from IZuckerman et all (Hp03), 
iKoester et alJ l)2005l) , or iBerger et all ll2005fl. Space ve - 
locities, v, relative to the sun are from |PauTf l|2003j) . 
IZuckerman et all l)2003[) . In eight cases where the space 
velocity is not known, we have replaced it with V3 
times the radial velocity, taken from the ci ted re f erence . 
Column d is the distan ce in pc, taken from IPaiilil l|2003l) . 
iMcCook fc SforlljlflP,^ . 



Object 



T cB logg [Ca/H] 





[K] 










[km/s] 


[pc] 


WD0032-175 


9235 


8 





-10 


5 


74 


32 


HS0047+1903 


16600 


7 


8 


-6 


1 


5 




HE0106-3253 


15700 


8 





-5 


8 


20 


66 


WD0208+296 


7201 


7 


9 


-8 


8 


89 




WD0235+064 


11420 


7 


9 


-9 


2 


11 


34 


WD0243-026 


6798 


8 


2 


-9 


8 


52 


21 


WD0245+541 


5190 


8 


2 


-11 


7 


38 


10 


HS0307+0746 


10200 


8 


1 


-7 


1 


38 




WD0408-041 


14400 


7 


8 


-6 


6 


40 


74 


WD0543+579 


8142 


8 





-10 


5 


120 


30 


WD0846+346 


7373 


8 





-9 


4 


20 


27 


WD1015+161 


19300 


7 


9 


-5 


9 


47 


95 


WD1102-183 


8026 


8 





-10 


3 


45 


40 


WD1116+026 


12200 


7 


9 


-6 


5 


26 


43 


WD1124-293 


9700 


8 


1 


-8 


2 


65 


34 


WD1150-153 


12800 


7 


8 


-6 





14 


85 


WD1202-232 


8800 


8 


2 


-9 


8 


37 


10 


WD1204-136 


11200 


8 





-7 


2 


56 


52 


WD1208+576 


5830 


7 


9 


-10 


8 


65 


20 


HE1225+0038 


9400 


8 


1 


-9 


7 


38 


29 


WD1257+278 


8481 


7 


9 


-8 


1 


50 


35 


HE1315-1105 


9400 


8 


4 


-9 


4 


17 


40 


WD1337+705 


20435 


7 


9 


-6 


7 


53 


33 


WD1344+106 


6945 


8 





-11 


3 


79 


18 


WD1407+425 


9856 


8 





-9 


8 


21 


26 


WD1455+298 


7366 


7 


6 


-9 


3 


101 


35 


WD1457-086 


20400 


8 





-6 


1 


24 


117 


WD1614+160 


17400 


7 


8 


-7 


2 


44 


117 


WD1633+433 


6569 


8 


1 


-8 


6 


30 


15 


WD1821-131 


7029 


8 


4 


-10 


9 


89 


20 


WD1826-045 


9200 


8 


1 


-8 


6 


26 


25 


WD1858+393 


9470 


8 





-7 


8 


3 


35 


WD21 15-560 


9700 


8 


1 


-7 


6 


54 


22 


HS2132+0941 


13200 


7 


7 


-7 


1 


37 


79 


WD2149+021 


17300 


7 


9 


-7 


7 


4 


25 


HE2221-1630 


10100 


8 


2 


-7 


2 


25 


52 


HS2229+2335 


18600 


7 


9 


-5 


9 


35 


91 


WD2326+049 


12100 


7 


9 


-6 


4 


59 


19 



— the H-rich white dwarfs observed with metals are cur- 
rently accreting in a steady state, with accretion and 
diffusion proceeding at the same rate; 

— there are no large-scale differences (like the two phases 
of the standard model). Many, and perhaps all, white 
dwarfs in the Local Bubble are accreting with a con- 
tinuous distribution of accretion rates caused by fluc- 
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Fig. 2. Diffusion timescales (in years) of Ca. The solid line 
indicates log g= 7.5, then following log g— 8.0, 8.5, 9.0 



tuations in ISM densities and differences in the stellar 
parameters. The ones we can identify as DAZs are the 
top of the distribution. 

Obviously this and the Dupuis et al. hypothesis are 
simplifications at extreme ends of the possibilities and 
could be characterized as "two-phase" vs. "continuous" 
accretion. In the following sections we explore whether 
the observations can lead to a distinction between the al- 
ternatives. 



3. The observed sample 

The data come from two large surveys: a search for Ca 
in DAs using the Keck telescope and the ESO Supernova 
la Progenitor Survey (SPY). Th e instruments and reduc- 
tion procedures ar e descr ibed in IZuckerman et "all l|2003h 
and IKoester et, all l(2005h . Atmospheric parameters T e ff, 
log g , and Ca abunda nces are given in those papers and 
in lBerger et al . (2005). Taking both sources together, we 
have a sample of 38 DAZs (excluding those marked as 
doub le degenerates or DA+DM pairs in lZuckerman et all 
2003) with observed photospheric Call lines and about 
450s DA with upper limits for the Ca abundance. 

Table [U collects the data for the DAZs, while Fig. [l] 
graphically shows the distribution of Ca abundances with 
T e g. The line corresponds to an equivalent width of the 
Call K line of 15 mA, which is the observational limit for 
our spectra with the highest signal-to-noise ratio. Because 
of the ionization of Call to Calll with increasing effective 
temperature, the lower limit for observable Ca increases 
until it becomes unobservable between 20000 and 25000 K. 
For all objects with Ca lines in our sample, the photo- 
spheric nature (as opposed to interstellar absorption) was 
confirmed through a comparison of radial velocity deter- 
minations from the Ca and the hydrogen Balmer lines. 

Ca is found at all temperatures filling the area of abun- 
dances between the lower observational limit and some up- 
per limit [Ca/H]« —5.7 ([Ca/H] is the logarithmic number 
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ratio of the elements) . This upper limit decreases towards 
lower T e g, which will be discussed in a later section. 

4. Diffusion timescales 

When trying to interpret the observational results in terms 
of an accretion-diffusion scenario, the most important the- 
oretical ingredient is the diffusion timescale. In the ab- 
sence of an outer supply, heavy elements disappear expo- 
nentially out of the reservoir (either the atmosphere or the 
outer convection zone, which is assumed homogeneously 
mixed) on a timescale given by (see e.g. iDunuis et all 
1992) 



TD 



M W p q 

Airr^pw 



(1) 



where My/D is the mass of the white dwarf, q the mass 
fraction of the reservoir, the radius at the bottom of the 
reservoir (always close to the white dwarf radius), p the 
mass density, and w the relative diffusion velocity between 
the heavy trace element and the main constituen t, which 
is hydrogen in our case. IPaauette et all l|l986b|) discuss 
the influence of the parameters q, p, w on the timescales, 
in particular with regard to the large differences between 
He and H envelopes. Basically the reason is the following: 
the convection zones (q) are much thicker in He envelopes 
at the same mass and effective temperature. The density 
at the bottom increases, but not as fast as q, and in addi- 
tion the velocity w is smaller in denser environments. The 
combination of these factors adds up to give timescales for 
He several orders of magnitudes larger than for H at the 
same T e g. 

The Dupuis et al. papers only give timescales in helium 
atmospheres in a tabular form, and only a few calculations 
for hydrogen envelopes below 10000 K, d ue to the lack of 
observa tional data at that time. Although IPaauette et alJ 
1986bh give data for hydrogen envelopes up to 15000 K, 



we nevertheless decided to repeat such calculations in or- 
der to obtain a homogeneous set of data covering the 
range of observations and incorporating the improved un- 
derstanding of the mixing-length parameters in DA white 
dwarfs. 

For these calculations we used our stellar atmosphere 
and stellar structure codes, which use up-to-date input 
physics for opacities and non-ideal gas effects in the equa- 
tion of state. The most relevant fact for this study is the 
use of the mixing- length version ML2 instead of ML3, 
and a mixing length of 0.6 pressure scale heights. Various 
mixing-length approximations used in white dwarf model- 
ing differ i n the choice of three numerical constants of the 
model fsee lFontaine et allll98lHTassoul et alTll990L for a 
definition of the nomenclature). In short, the ML3 version 
is a more efficient convection with higher energy flux at the 
same temperature gradient. Necessary stellar parameters 
for the diffusion time calculations, e.g. the physical condi- 
tions at the bottom of the reservoir of heavy elements to 
be depleted, were obtained from these models. Convection 
zones start very shallow in DAs above T e g« 15000 K, and 



the reservoir is the atmosphere. We define, somewhat ar- 
bitrarily, the limit at an optical depth of tr oss = 5. When 
a 0.6 Mq white dwarf cools down below 13000 K, the bot- 
tom becomes deeper than this, and the reservoir is defined 
by the bottom of the convection zone. For a more (less) 
massive white dwarf, the transition occurs at a slightly 
higher (lower) T e g. 

Collision integrals fo r the diffusion v e locities were ob- 
tained from the fits in IPaauette et al.1 I \ 986al) and the 
calcula tions followed the prescriptions in 



Paauette et alJ 



(1986b) very closely. One major exception was that the 
detailed distribution over ionization states was obtained 
from the models and average diffusion velocities deter- 
mined instead of just using the dominant ions as in the 
latter work. In view of significant uncertainties in non- 
ideal gas effects (e.g. pressure ionization) or convection 
theory, we do not believe that the result will be more ac- 
curate, but this method avoids discontinuities in timescale 
changes with temperature of the models. Table gives 
the resulting timescales for Ca, Mg, and Fe in hydro- 
gen, while Fig. El shows the Ca data in graphic form. In 
ge neral these resu l ts are very similar to tho se obtained 
bv IPaauette et all l|l986b|) and earlier work l|MuchmoreL 
Il979lll984l) .' The main difference is that the steep increase 
in diffusion timescales towards cooler temperatures occurs 
at slightly lower T c g because of different formulations of 
the mixing-length theo ry used, e.g. ML3 in the case of 
IPaauette et al.l l|l986b|) . 

The first conclusion from these data is a confirmation 
of earlier results: the diffusion timescales are extremely 
short compared to evolutionary timescales, and also much 
shorter than diffusion timescales in He-rich envelopes. 
Over most of the observed range, they are much shorter 
than 1000 yr, going down to days at the high-temperature 
end; and only for the very coolest objects do the timescales 
approach 100000 yr. The steep increase in the timescales 
between 12500 and 10000 K is directly connected to the 
very strong increase in the depth of the convection zone, 
when the white dwarf cools through this temperature re- 
gion. If the distribution in Fig.^is due to objects after the 
end of their accretion episode (note that an exponential 
decline of a metal abundance means equal time for each 
decade of abundances and thus a "homogeneous" popula- 
tion of the area above the visibility limit), then the metals 
will be visible only for a very short time compared to the 
cooling time through this range - approximately 5.6 x 10 s 
yr from 200 00 to 10000 K, and 1.7 x 10 9 yr from 10000 
to 6000 K l(Woodl . Il995h . This conflicts with the result 
that a significant fractio n of DA shows me tals, > 20% in 
IZuckerman et alJ l)2003[) and still > 5% in iKoester et alJ 
(2005), a sample much more biased towards higher T e g. 
This high fraction is only possible with a constant supply 
of new DAZs through accretion. On the other hand, tak- 
ing the numbers for the two-phase scenario at face value, 
where the star spends on the order of 10 6 yr in the accre- 
tion phase, one would expect a completely different dis- 
tribution with the vast majority of the objects near the 
upper boundary of abundances. This leaves little doubt 
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Fig. 3. Accretion rates for 38 DAZ in solar masses per 
year, assuming that the Ca abundance in the accreted 
material is solar. The line shows the limit of the accretion 
rates resulting from an observational limit of 15mA for a 
0.6 M© white dwarf 

that this interpretation is not correct and that we have to 
assume that accretion is still ongoing in most, if not all 
objects. This conclusion is inevitable for the hot objects 
with timescales of a few days (e.g. about 10 days for G29- 
38, which has been observed to have Ca since 1997), but 
is very plausible for all objects. We will thus assume that 
all DAZ white dwarfs are currently accreting at a steady 
state (between accretion and diffusion), with this gener- 
aliz ation made plausible by the short diffusion timescales 
(see lDutmis et al 1 11 9921 for a discussion). 

5. Steady state accretion rates 

With the assumption of steady state acc retion, we can 
immed iately determine the accretion rate (jPupuis et all 

M acc = M^R^L (2) 

TD-^acc 



Table 3. Diffusion timescales td, mass fraction at the 
layer of diffusion, logg, and the accretion rates M in 
M /yr for 38 DAZs abundances. 



Object log-To Mwd log? logM 



i ti~i r\ r\ n r\ n "7 rr 

WD0032-175 


2 


72 





600 


-10 


04 


-17 


84 


HS0047+1903 


-1 


48 





512 


-15 


56 


-14 


83 


HE0106-3253 


-1 


82 





616 


-15 


90 


-14 


45 


WD0208+296 


3 


78 





537 


-8 


43 


-15 


64 


unmoor i r\f* a 

WD0235+064 


-1 


04 





550 


-14 


91 


-17 


69 


I TT"t AO/ID AAC 

WD0243-026 


3 


39 





719 


-8 


67 


-16 


37 


WD0245+541 


4 


47 





712 


-6 


98 


-17 


66 


TTO AOA7 i f\T A f 

HS0307+0746 


1 


81 





664 


-11 


42 


-14 


87 


I Tl~\ A/1 AO A/1H 

WD0408-041 


-1 


56 





506 


-15 


70 


-15 


39 


T TT^i AT /I A i I - *7/"\ 

WD0543+579 


3 


22 





597 


-9 


17 


-17 


48 


I TI~t AO/1C i A /I A 

WD0846+346 


3 


51 





595 


-8 


72 


-16 


22 


I TT"t A f\ A FT i A H 

WD1015+lbl 


-1 


59 





571 


-15 


60 


-14 


51 


WD1102-183 


3 


27 





597 


-9 


10 


-17 


25 


WD1116+026 


-1 


73 





552 


-15 


86 


-15 


24 


WD1124-293 


2 


29 





663 


-10 


66 


-15 


69 


WD1150-153 


-1 


63 





500 


-15 


79 


-14 


82 


WD1202-232 


2 


59 





724 


-9.900 


-16 


79 


T TT^ H A A /I A O 

WD1204-136 


-0 


34 





606 


-14 


01 


-15 


44 


WD120o+576 


4 


50 





530 


-7 


31 


-17 


25 


HE1225+0038 


2 


48 





662 


-10 


31 


-17 


03 


T TT^ A AT7 i O ^7 O 

WD1257+278 


3 


26 





541 


-9 


25 


-15 


24 


T TT 1 -lOHCr H 4 AT 

HE1315-1105 


2 


02 





857 


-10 


63 


-16 


47 


I Tl~t A D n "7 i 7 A r 

WD1337+705 


-1 


57 





574 


-15 


56 


-15 


29 


WD1344+106 


3 


68 





594 


-8 


46 


-18 


02 


WD1407+425 


2 


27 





602 


-10 


83 


-17 


48 


WD1455+298 


4 


33 





386 


-7 


98 


-16 


38 


WD1457-086 


-1 


73 





628 


-15 


71 


-14 


64 


I Tl~t A f A A i A C r\ 

WD1614+1dO 


-1 


46 





514 


-15 


53 


-15 


91 


Wi/ J. doo~*±oo 


3 


69 


o 


654 


-8 


34 


-15 


17 


WD1821-131 


2 


89 





853 


-9 


21 


-17 


43 


WD1826-045 


2 


58 





661 


-10 


11 


-15 


83 


WD1858+393 


2 


58 





601 


-10 


30 


-15 


27 


WD2115-560 


2 


29 





663 


-10 


66 


-15 


09 


HS2132+0941 


-1 


47 





454 


-15 


63 


-15 


97 


WD2149+021 


-1 


62 





565 


-15 


68 


-16 


37 


HE2221-1630 


1 


84 





727 


-11 


28 


-14 


82 


HS2229+2335 


-1 


60 





569 


-15 


62 


-14 


53 


WD2326+049 


-1 


70 





551 


-15 


81 


-15 


13 



with mass fraction of the heavy element in the reservoir 
X and its abundance in the accreting matter X acc . We 
assumed the solar abundance of Ca for X acc and took the 
mass fraction q of the reservoir from our models with the 
appropriate stellar parameters. The resulting mass accre- 
tion rates are given in Table Eland Fig. Tue radius does not change significantly, but the increased 
The continuous distribution of observed Ca abun- diffusion flux pw leads to lower abundances in the ob- 
dances thus translates into a continuum of derived ac- servable reservoir. The range o f accretion rat e s is in g ood 
cretion rates, with a maximum rate of approximately agreement with the results of |Dupuis e£aQ <|l993bf) for 
3 x 10~ 15 M Q /yr over the whole range of observed tern- helium-envelope white dwarfs; in fact, their "high accre- 
peratures. The lack of higher abundances below T eff w tion rat e" used for the comparison with observed abun- 
10000 K in Fig. [J is due to the bottom of the reservoir dances is 5 x 10~ 15 M /yr, in excellent agreement with 
moving into deeper and denser regions, as can be seen Fig. El If the maximum accretion rates are indeed indc- 
from the combination of Eqs.(l) and (2): pendent of the surface temperature of the star, one could 

even speculate that our mixing-length description predicts 

M acc X acc Tu M acc X acc convection zone depths at the low temperature end, which 

~ M WD q ~ Anrlpw ' are slightly too large! 
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Table 2. Logarithm of the diffusion timescales (in years) for Ca, Mg, and Fe in hydrogen- rich white dwarfs as a 
function of effective temperature and surface gravity 



T cff [K] 




7.5 




7.75 




8.0 


logs 
8.25 




8.5 




8.75 




9.0 




Ca 


6 


1353 


5 


7436 


5 


3045 


4 


7835 


4 


1353 


3 


4683 


2 


5939 


5000 


Me 


6 


1629 


5 


7188 


5 


2500 


4 


8033 


4 


2051 


3 


5215 


2 


6560 




Fe 


6 


1406 


5 


7159 


5 


2216 


4 


6712 


3 


9877 


3 


3064 


2 


4714 




Ca 


5 


1778 


4 


6687 


4 


2179 


3 


8380 


3 


4399 


3 


0440 


2 


5471 


6000 


Me 


5 


2526 


4 


7300 


4 


2702 


3 


8767 


3 


4632 


3 


0424 


2 


5036 




Fe 


5 


1932 


4 


6841 


4 


2439 


3 


8780 


3 


4981 


3 


1134 


2 


5948 




Ca 


4 


6835 


4 


1794 


3 


6554 


3 


1968 


2 


7003 


2 


2280 


1 


7375 


7000 


Me 


4 


7280 


4 


2117 


3 


6735 


3 


2086 


2 


6994 


2 


2151 


1 


7150 




Fe 


4 


6304 


4 


1124 


3 


5628 


3 


0924 


2 


5650 


2 


0642 


1 


5538 




Ca 


4 


2452 


3 


7668 


3 


2781 


2 


8298 


2 


3413 


1 


8960 


1 


3673 


8000 


Me 


4 


2724 


3 


7838 


3 


2839 


2 


8307 


2 


3366 


1 


8874 


1 


3521 




Fe 


4 


1262 


3 


6333 


3 


1189 


2 


6570 


2 


1427 


1 


6812 


1 


1178 




Ca 


3 


6618 


3 


2656 


2 


8340 


2 


4193 


1 


9883 


1 


5547 


1 


0580 


9000 


Me 


3 


6729 


3 


2722 


2 


8328 


2 


4146 


1 


9811 


1 


5451 


1 


0462 




Fe 


3 


4729 


3 


0716 


2 


6270 


2 


2022 


1 


7571 


1 


3075 





7892 




Ca 


2 


1235 


2 


1596 


2 


1125 


1 


8695 


1 


5642 


1 
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7282 


10000 


Me 


2 


0589 


2 


1135 


2 


0940 


I 
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1 


5580 


1 


1841 
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Fe 
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9373 


1 


9567 


1 


9125 


1 


6605 


1 


3344 





9397 





4499 




Ca 





7272 


1 


0429 


1 


3574 


1 


3264 


1 


1629 





9366 





5481 


10500 


Me 





7107 





9840 


1 


2991 


1 


2851 


1 


1402 





9232 





5323 




Fe 





6606 





9191 


1 


1606 


1 


1139 





9485 





7065 





2875 




Ca 


-0 


6754 


-0 


2462 





2184 





4745 





5944 





5090 





2723 


11000 


Me 


-0 


4665 


-0 


1343 





1712 





4423 





5441 





4723 





2487 




Fe 


-0 


8266 


-0 


2525 





1212 


o 


3335 





3868 





2845 





0385 




Ca 


-1 


0396 


-1 


2166 


-1 


0994 


-0 


5974 


-0 


2729 


-0 


1556 


-0 


1807 


11500 


Me 


-0 


8548 


-1 


0199 


-0 


8862 


-0 


5751 


-0 


3171 


-0 


1855 


-0 


2240 




Fe 


-1 


2632 


-1 


3974 


-1 


1010 


-0 


6656 


-0 


4005 


-0 


3263 


-0 


4012 




Ca 


-1 


1486 


-1 


4876 


-1 


7497 


-1 


6973 


-1 


2673 


-0 


9217 


-0 


8417 


12000 


Me 


-0 


9675 


-1 


3054 


-1 


5580 


-1 


4801 


-1 


1788 


-0 


9544 


-0 


8612 




Fe 


-1 


3827 


-1 


7171 


-1 


9453 


-1 


7734 


-1 


3084 


-1 


0355 


-0 


9847 




Ca 


-1 


1597 


-1 


5499 


-1 


9446 


-2 


2830 


-2 


5642 


-2 


8837 


-2 


7075 


13000 


Me 


-0 


9673 


-1 


3627 


-1 


7644 


-2 


0971 


-2 


3796 


-2 


6888 


-2 


5069 




Fe 


-1 


3675 


-1 


7691 


-2 


1787 


-2 


5123 


-2 


7764 


-3 


0656 


-2 


7119 




Ca 


-1 


1156 


-1 


5065 


-1 


8966 


-2 


2858 


-2 


6730 


-3 


0616 


-3 


3721 


14000 


Me 


-0 


9063 


-1 


3008 


-1 


6961 


-2 


0920 


-2 


4880 


-2 


8880 


-3 


1820 




Fe 


-1 


2870 


-1 


6861 


-2 


0857 


-2 


4867 


-2 


8887 


-3 


2975 


-3 


5648 




Ca 


-1 


0700 


-1 


4602 


-1 


8491 


-2 


2354 


-2 


6268 


-3 


0241 


-3 


4227 


15000 


Me 


-0 


8583 


-1 


2519 


-1 


6451 


-2 


0368 


-2 


4286 


-2 


8253 


-3 


2286 




Fe 


-1 


2172 


-1 


6147 


-2 


0126 


-2 


4102 


-2 


8097 


-3 


2097 


-3 


6109 




Ca 


-1 


0351 


-1 


4223 


-1 


8116 


-2 


2034 


-2 


5981 


-2 


9955 


-3 


3929 


16000 


Me 


-0 


8205 


-1 


2101 


-1 


6031 


-1 


9952 


-2 


3871 


-2 


7870 


-3 


1872 




Fe 


-1 


1619 


-1 


5572 


-1 


9569 


-2 


3555 


-2 


7510 


-3 


1483 


-3 


5439 




Ca 


-1 


0070 


-1 


3937 


-1 


7838 


-2 


1772 


-2 


5744 


-2 


9721 


-3 


3693 


17000 


Me 


-0 


7924 


-1 


1822 


-1 


5741 


-1 


9612 


-2 


3587 


-2 


7580 


-3 


1572 




Fe 


-1 


1056 


-1 


5030 


-1 


9037 


-2 


2994 


-2 


6988 


-3 


0968 


-3 


4931 




Ca 


-0 


9846 


-1 


3722 


-1 


7646 


-2 


1596 


-2 


5554 


-2 


9530 


-3 


3503 


18000 


Me 


-0 


7697 


-1 


1600 


-1 


5492 


-1 


9401 


-2 


3369 


-2 


7358 


-3 


1344 




Fe 


-1 


0435 


-1 


4465 


-1 


8485 


-2 


2492 


-2 


6490 


-3 


0523 


-3 


4543 




Ca 


-0 


9681 


-1 


3565 


-1 


7488 


-2 


1438 


-2 


5397 


-2 


9373 


-3 


3346 


19000 


Mg 


-0 


7512 


-1 


1406 


-1 


5288 


-1 


9234 


-2 


3199 


-2 


7184 


-3 


1167 




Fe 


-0 


9809 


-1 


3866 


-1 


7916 


-2 


1955 


-2 


6010 


-3 


0076 


-3 


4125 




Ca 


-0 


9545 


-1 


3433 


-1 


7352 


-2 


1304 


-2 


5263 


-2 


9239 


-3 


3213 


20000 


Mg 
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-1 
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-1 


5148 


-1 
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-2 
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-2 
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Fe 
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-1 
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-2 


5525 


-2 
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-3 


3692 




Ca 
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8771 


-1 
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-1 


6684 


-2 


0677 


-2 
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-2 


8662 


-3 


2657 


24000 


Mg 


-0 


7030 


-1 


0903 


-1 


4830 


-1 


8772 


-2 


2719 


-2 


6675 


-3 


0640 




Fe 


-0 


7958 


-1 


1973 


-1 


6004 


-2 


0050 


-2 


4111 


-2 


8187 


-3 


2277 
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6. Correlation of accretion rates with stellar 
parameters 

What is the reason for the large variation in accretion 
rates, e.g. by three orders of magnitude at the cool 
end? To study this question, we have to proceed to 
the theoretically much less secure ground of the accre- 
tion process. Rather than discuss the problems of the 
physics of accretion, which is beyond the scope of the 
present study, we assume that the mass accretion rate is 



given by the Bondi-Ho yle formulationj|Ho i vjB i ^ i I i vttle^or. 
1933 iBoudi fc HovleL 17941 iRondl Il9fi* iRnffertl llflfr 



Ruffert fc ArnettJ . [l994) : 



M ac 



( c 2 +u 2)3/2 



(4) 



Here, M is the stellar mass, v s its velocity relative to the 
interstellar gas at large distance, c the sound velocity and 
p the unperturbed ISM density. For the warm (ss 8000 K) 
phase of the ISM, the sound velocity is not important com- 
pared to the space velocity for most objects; we have taken 
c = 8 km/s and use a new variable v = y^c 2 + v 2 in the 
following. The accretion rate therefore can be separated 
into individual properties of the stars (oc M 2 /v 3 ) and the 
condition of the ambient medium given by p. Fortunately, 
the stellar parameters mass and space velocity (in 8 ob- 
jects estimated from the radial velocity) are known for our 
complete sample. 

Figure 0] shows accretion rates as a function of the 
combination M 2 /v 3 ("stellar factor"), which enters the 
Bondi-Hoyle formula Eq. (4). Although the stellar factor 
varies by more than three orders of magnitude, there is 
obviously no correlation with the accretion rates. At most 
one could state that there is a small tendency towards 
higher rates for large factors and vice versa. This can only 
mean that the second factor, the interstellar density, must 
play an equally large role and show similar variations if 
accretion is coming from the ISM. This is demonstrated in 
Fig. [3 which shows the distribution of interstellar density 
(derived from all the other known factors in Eq. (4) and 
converted to hydrogen particle density) as a histogram. As 
expected from the lack of correlation in Fig. [3 the derived 
interstellar densities span the range rig = 0.01 — 1 cm~ 3 , 
with very few objects outside (we use nn for the sum of 
neutral and ionized hydrogen atoms). It should be noted 
that this number is derived from the mass density and 
thus includes neutral and ionized hydrogen. The direct 
(not log) average of all values is ss 0.5, which is very likely 
biased towards higher values, because for higher T e g we 
can observe only the highest accretion rates. 

For most of the current sample the distances are 
known, and we can therefore determine the location in the 
solar neighborhood corresponding to the individual den- 
sity determinations. We can also derive upper limits for 
the ISM densities for many of the numerous objects with 
upper limits for the Ca abundance (Fig. [I}. We had hoped 
to find examples of several DAZs close together in an area 
of the sky, all with enhanced ISM densities derived from 
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Fig. 4. Accretion rates vs. the "stellar factor" M 2 /v 3 for 
38 objects. Objects with T c ff< 10000 K are marked with 
circles. Accretion rates are in Mq per year, the units for 
the stellar factor are M Q and km/s. 
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Fig. 5. ISM density (neutral + ionized hydrogen per cm 3 ) 
derived from assuming the validity of the Bondi-Hoyle ac- 
cretion formula. See text 

the accretion rates, but no such correlations are apparent. 
If our interpretation up to this point is at all meaningful, 
this means that the ISM density varies by two orders of 
magnitude on the scale of a few pc. 

7. Discussion and conclusions 

The two-phase accretion scenario in the highly idealized 
form of the Dupuis et al. work cannot explain the ob- 
servations of the DAZs discussed in this study. The ab- 
sence of cold neutral matter in the solar neighborhood, 
the large fraction of DA showing metals in spite of short 
diffusion timescales, and the distribution of the observed 
abundances all point to ongoing accretion even under the 
conditions of the ISM in the Local Bubble. 

Although we show that the stellar parameters entering 
the calculation of accretion rates according to Bondi-Hoyle 
hydro dynamic accretion vary over several orders of mag- 
nitude, this alone cannot explain the variation in observed 
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accretion rates. Additionally, we need a significant varia- 
tion in the ISM densities, with typical values of ionized 
plus neutral hydrogen nji sa 0.1 — 0.2 cm -3 and variation 
by a factor of 10 in both directions. The only viable candi- 
date is the warm, partially ionized phase of the solar neigh- 
borhood ISM , which has b e en stu died recently in a series 
of papers bv lLinskv et all <j2000|) and Rcd field fc Linskvl 
fc00al200al2004albb . 

From these papers the following picture emerges. The 
absorption lines along the line-of-sight to many nearby 
stars show the existence of distinct components, ten- 
tatively called "clouds". Neutral hydrogen densities are 
typically tlhi ~ 0.1 cm" 3 , with hydrogen ions adding 
rip w 0.11 cm" 3 to the mass. The length scale is 2.2 pc, 
with a range of 0.1 to 11 pc. Our own sun lies within a 
similar cloud, the Local Interstellar Cloud or LIC, with ex- 
tensions of 4 - 6 pc. We use these numbers for some very 
simple estimates, to test the plausibility of a connection 
between these clouds and the DAZ phenomenon. 

Assuming current accretion within such a cloud to be 
responsible for the observed metals, the observations de- 
mand a filling factor on the order of 10% within 50 pc. 
Taking spherical clouds with a radius of 1.5 pc, the total 
number of clouds in this volume would be 3700. The total 
solid angle covered by these clouds (distributed homoge- 
neously) is 31.4 or 2.5 times the total sphere of 4ir. This 
is in excellent agreement with the findings of 1-3 com- 
ponents (and an average of w 2) along each line-of-sight 
in the Redfield and Linsky studies. A white dwarf with 
a space velocity of 30 km/s would travel 1 pc in about 
33000 yr, plenty of time to reach steady state abundances 
and stay constant over the observational timescales. 

The typical distance between the centers of any two 
clouds would be 3.2 pc, about the same as the size of 
the clouds. This should serve as a caveat and a reminder 
that the concept of "clouds" may again be an oversimpli- 
fication. However, this does not affect our interpretation, 
since all we need is a scale length for significant variation 
in the physical parameters of the ISM. These estimates are 
all quite favorable to the hypothesis of ongoing accretion 
from the warm phase of the ISM, alleviating the strongest 
objection against the accretion/diffusion scenario - the 
absence of cool, neutral clouds. 

There are, however, other open questions remaining. 
One of them is whether the assum ption of f l uid (B ondi- 
Hoyle) accretion rates is justified. iKoesterl l)l97fil) con- 
cluded that accretion rates would be much lower than the 
fluid rate since interactions are not sufficient to destroy 
the momentum perpendicular to the accretion column, a 
necessary condit ion for the hydrodynamic treatment. On 
the other hand. lAlcock fc Illarionovl 1 19801) argued that an 
ionized plasma would always accrete at the fluid rate due 
to plasma instabilities and magnetic fields. 

Another difficult question is the composition of the ac- 
creted matter. We have derived a very robust result on the 
accretion of Ca atoms (which can be derived from the mass 
accretion rates in Table by multiplying those numbers 
by 8 x 10 -5 ). The total accretion rate follows from as- 



suming a solar composition. There is no evidence for that; 
on the contrary, observations in helium-rich white dwarfs 
with metals indicate that in many cases very little or no 
hydrogen is accreted. This seems to show that the accre- 
tion is preferentially of dust grains, either interstellar or 
from some circumstellar cloud or disk. Circumstell ar ma- 
terial has been found for G29-38 and GD362 (see e.g.lJura , 
|2003 f IReach et all 120051: IBecklin et all 120051: iKilic et al. . 
2005, for recent studies), with the origin attributed to 
comets, tidal disruption of pla nets or asteroids, or the 
merging of two white dwarfs ijLivio et all l2005|) . It is 
quite possible that in such cases accretion from a different 
source than the ISM is responsible for the observed met- 
als; similarly, in several c ases of close M dwarf companions 
((Zuckerman et all l2003h , a stellar wind may provide the 
heavy elements. 

However, these seem to be excepti onal cases. As found 
in the comprehensive recent work bv iFarihi et alJ l|2005h 
none of the 371 white dwarfs studied showed a near- 
infrared excess indicative of circumstellar dust similar to 
G29-38 and GD362. One should notice that the search was 
limited by the signal-to-noise ratio of the Ks magnitude 
in the 2MASS sample and only about 1/3 of the sample 
had adequate data, so further detections may be forth- 
coming from ongoing searches (J. Farihi, priv. comm.). 
Nevertheless, it appears unlikely that such processes can 
account for the 10 - 20% of hydrogen-rich white dwarfs 
that are below 20000 K and have heavy elements. 
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